Periodic Orbits and Binary Collisions in the Classical Coulomb Three-Body Problem 
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In the helium case of the classical Coulomb three-body problem in two dimensions with zero 
angular momentum, we develop a procedure to find periodic orbits applying two symbolic dynamics 
for one-dimensional and planar problems. A sequence of periodic orbits are predicted and are 
actually found numerically The results obtained here will be a cornerstone for finding the remaining 
periodic orbits, which needed for semiclassical applications such as periodic orbit quantization. 
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The microscopic three-body problem, for example, the 
helium atom, attracts both theoreticians and experimen- 
talists because of rich structure of its spectrum. In par- 
ticular, the excited states below the double ionization 
threshold (E = 0) show the most complicated spectral 
structure, and contain various information on electron- 
electron correlations. This spectral region is not only 
most interesting but also most difficult in both theoret- 
ical and experimental aspects. In investigating its spec- 
trum, semiclassical methods inevitably face with chaotic 
nature of its classical dynamics pQ. In a theoretical as- 
pect, for the collinear eZe configuration, the semiclassi- 
cal periodic orbit quantization was carried out by us- 
ing the celebrated Gutzwiller trace formula [2] within 
reasonable accuracy [3] [4]. However, such successes 
of semiclassical and classical approaches have been re- 
stricted to one-dimensional case, namely the collinear eZe 
and eeZ configurations and the Wannier ridge configura- 
tion [TJ |3j HI [5]. The helium atom in two dimensions 
remains as a challenging object. 

Main difficulties are due to (1) the mixture of chaos 
and tori, (2) the high-dimensionality, and (3) the sin- 
gular nature of its dynamics. To overcome these dif- 
ficulties, one needs sophisticated mathematical tools. 
McGehee's blow-up transformation and the concept of 
the triple collision manifold [6] are important exam- 
ples H U EB HB CE2 CES]. Thanks to these tools, the 
dynamics near triple collision and the structure of sta- 
ble and unstable manifolds were elucidated to some ex- 
tent. In fact, a quantum manifestation of triple collision 
was demonstrated [14j [15] . But we still know less about 
periodic orbits in two-dimensional dynamics, which are 
needed for semiclassical analysis. 

Quite recently, one additional tool, i.e., the symbolic 
dynamics for the planar gravitational problem is devel- 
oped [16 . In this method, binary collisions are naturally 
obtained as boundaries of different symbol sequences. 
We already have one-dimensional symbolic dynamics [8] . 
Then, we can systematically search periodic orbits in 



the dynamics with zero initial velocities (DZIV) using 
one-dimensional and two-dimensional symbolic dynam- 
ics. We do this in the present report. 

As is known in [10], the initial condition space of the 
DZIV is bounded by the collinear initial conditions (see 
below). First, using the two-dimensional symbolic dy- 
namics, we numerically confirm the existence of curves 
of initial conditions of orbits exhibiting binary collisions. 
We call these the binary-collision curves (BCCs). Next, 
we go back to symbolic dynamics for the collinear prob- 
lem and obtain symbol sequences of triple collision orbits. 
As predicted in [10], a sequence of BCCs is shown to have 
two end-points on a = 7r. In other words, these BCCs 
can be uniquely specified by the end-points. Then, we 
informally consider a one-dimensional symbolic dynam- 
ics along BCCs and obtain self-retracing periodic orbits 
along them. Our procedure (Fig. [3| will be justified by 
numerical integration (Fig. [4|. Finally, we characterize 
periodic orbits with physical quantities. 

Now we introduce the system considered. In the he- 
lium atom, two electrons are denoted by particles 1 and 
2, and the nucleus is denoted by particle 3. Let be the 
position vector of the ith particle, and let = ri — rj. 
We assume that the nucleus has infinite mass and that 
initial velocities of particles are zero, which implies that 
the problem is planar. Then, the Hamiltonian in the 
hyperspherical coordinates is given by 

H = 1 (r? I I \ 

2\ Pr ^ r 2 + r 2 sin 2 ( X ) J 



r ycos(|) sin (|) Z[l — sin(x) cos(a)] 1 / 2 

where r = (r 2 + r^) 1 ^ 2 is the hyperradius with r\ — 
rcos(x/2), T2 = rsin(x/2), and x = 2arctan(r2/ri). 
Z(= 2 in this Letter) is the charge of the nucleus. Z ap- 
pears in the denominator of the potential term because 
of the scaling of variables, a is the angle between the 
vectors ris and r^z- The total energy is E = H. 
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Next we introduce symbolic dynamics [16] . Three par- 
ticles from a triangle. The area of this triangle has the 
sign r\ = e z • (ri3 x r 2 3)/|^i3 x ^23 1- The length of the 
ith side of the triangle is ^ = \vjk\, where k) is a 
cyclic permutation of (1,2,3). If 77 changes its sign from 
plus to minus and max{^,/2,fe} = h, then we assign a 
symbol k. If rj changes its sign from minus to plus and 
max{/i, Z2, ^3} = h, then we assign a symbol fc + 3. Then 
we obtain a symbol sequence s = •S1S2S3 * * ' for a given 
orbit with the symbol set S = {1, 2, 3, 4, 5, 6}. 

By scaling the system, we remove the dependence on r 
from the DZIV. Then the initial conditions are uniquely 
specified by < x < n an d < a < 2tt. Moreover, the 
quarter of the initial condition space D 1 / 4 = {(%, a)|0 < 
X^f^O^^^ 71 "} suffices for our purpose owing to 
the symmetry of the problem. We integrate orbits start- 
ing at points in .D1/4, and assign a symbol sequence to 
each orbit following the procedure given in [16]. Then, 
is partitioned into regions of different symbol se- 
quences. One example of partitions by symbol sequences 
up to symbol length 8 is shown in Fig. [I] The BCCs are 
clearly seen as the boundaries of the partitioned regions. 
This shows a big advance compared with [10 . In that 
work, BCCs were obtained point by point along them, so 
the task was rather time consuming. In addition, orbit 
integrations near a = tt were difficult due to the closeness 
to triple collision. The present method basically avoids 
triple collision, so is accurate and fast. 

The intersections of BCCs with a = ir are the ini- 
tial conditions for triple-collision orbits (TCOs). We call 
them the triple-collision points (TCPs). We attach, in 
Fig. [l] the orbits of some of the TCPs with symbol se- 
quences 0, 210, 20, and 220 where '0' represents triple col- 
lision, '1' binary collision between particles 1 and 3, and 
'2' binary collision between particles 2 and 3. Hereafter 
the symbolic dynamics means the symbolic dynamics in 
the collinear eZe configuration. 

We observe two groups of BCCs. One group, which 
we call the first kind, consists of the BCCs which start 
at a = 7r and end at a = n. We denote a BCC of this 
group by a lobe from its shape, and express it as w m 0- 
w n using the symbol sequences w m and w n of the left 
and right end points where ^ is a word of length k. 
The other, which we call the second kind, consists of the 
BCCs which start at a = tt and end at a = 0. There 
are other possibilities of the terminals of the BCCs. For 
our present purpose, we only consider periodic orbits on 
BBCs of the first kind. There exist an infinite number of 
TCPs on a — tt. Correspondingly, the number of lobes 
is infinite (see also Fig. [5] and [8 ), which immediately 
implies that there exist infinitely many (unstable) POs 
on BCCs of the first kind, as we show later. 

The appearance order of TCPs is illustrated in Fig. [2] 
If the symbol length is the total number of the TCPs 
appearing on a = tt is 2 l ~ x where / denotes the gener- 
ation of TCPs. The successive generations of TCPs are 
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FIG. 1: Partition in D1/4 and triple-collision points in 
the collinear eZe configuration: The large figure represents 
the partition obtained from the symbol sequences of symbol 
length 8 in the symbolic dynamics for two-dimensional dy- 
namics. The boundaries between the regions of the partition 
are binary-collision curves. The small figures represent the 
triple-collision orbits in the collinear eZe configuration, (a) 0. 
(b) 210. (c) 20. (d) 220. The blank region at the lower-right 
corner is the forbidden region (E > 0). 
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FIG. 2: Appearance of triple-collision points on the line 
a — 71 : I is the symbol length in the symbolic dynamics for 
the collinear eZe dynamics. The symbol sequences of triple- 
collision orbits on a = ty are represented in the symbolic dy- 
namics of the collinear eZe case. The solid (resp. dotted) line 
represents the lobe corresponding to 2-3 (resp. 1-3) collision. 



connected by the following rule: In the (I + l)th genera- 
tion, new TCPs are added in between the already existing 
TCPs in the Ith generation. Let w be a word of length 
I — 1, and wO be a symbol sequence of a TCP which al- 
ready exists in the Ith generation. Then, wlO appears 
to the immediate right of w0, while w20 appears to the 
immediate left of wO. (The former rule does not apply 
for I = 1.) Any orbit in BCCs with this new wlO or w20 
at the right (resp. left) end experiences a 2-3 (resp. 1-3) 
collision. 

Here we give an intuitive proof of the existence of a PO 
on each BCC of the first kind (i.e., lobe). We introduce 
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FIG. 3: Tails of orbits on the lobe w20-w210 



coordinate ft (0 < ft < 1) on the curve so that ft = 
at the left end and ft = 1 at the right end. The orbit 
changes its topology when we move along the curve from 
ft = to ft = 1 since the orbits at both ends have different 
histories of collisions. There are three types of lobes: (i) 
lobe w20-w0 (w ^ 2 k , k > 0) or lobe wO-wlO; (ii) lobe 
wO-w'O where symbol lengths of w and w' differ by 2 or 
more; (hi) lobe 2 /e+1 0-2 /c (k > 0). 

First consider case (i), i.e., the case of lobe w20-w210. 
Let us move from w20 (ft = 0) to w210 (ft = 1). The rel- 
ative position of the minimum of r\ continuously changes 
as in Fig. [3j The triple collision at ft = becomes a 2-3 
binary collision for < ft < 1. As ft increases, the next 
minimum of T\ approaches the 2-3 binary collision from 
the right, and finally coincides with it at ft = 1 resulting 
in triple collision. In this process, there exists a parame- 
ter value such that the 2-3 binary collision and the local 
maximum of r\ take place at the same time (see the third 
panel of Fig.|3|. At this moment, the angular momentum 
of electron 2 and the nucleus is zero with respect to their 
center of mass, which means that the angular momentum 
of electron 2 is zero, since the nucleus has infinite mass. 
This in turn means that the angular momentum of elec- 
tron 1 is zero, which implies that electron 1 stands still 
at this moment. Then, both electrons retrace the path 
they tread, that is, a self-retracing PO is obtained. In 
the above proof, the persistence of the minima of T\ and 
V2 is crucial. 

The proof for case (ii) may be similar to case (i). How- 
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FIG. 4: (Color online) Periodic orbits (POs) found : The 
nucleus is fixed at the origin. For the POs with binary colli- 
sions, names of the lobe to which they belong are indicated. 
In the orbit integration, we set E = — 1. 



ever, there is a possibility that the orbit escapes, that is, 
one of electrons escapes to infinity. We skip this case, 
since we need a special care to treat the non-persistence 
of minima of T\ and/or T2- 

In case (iii) , the situation is different from that of case 
(i). In this case, we can draw a similar orbital change as 
in Fig. [3] This time, however, the third minimum of T2 
disappears due to the escape of the orbit when the first 
minimum approaches the second minimum of T2- Then 
we do not have the simultaneous occurrence of the max- 
imum of r2 and the minimum of 7*1 . Therefore, there is 
no POs on the 2 /c+1 0-2 fe lobes. 

To obtain POs numerically, we look for the moment 
of zero velocities other than the initial moment. We ac- 
tually carried out this procedure with many trials and 
errors. Then we find some POs using the improved New- 
ton method. 

We show example trajectories of POs in Fig. |1J We 
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FIG. 5: The Lyapunov exponents and the period for a family 
of periodic orbits: (a) The Lyapunov exponents, (b) The 
period. n r is the number of revolutions. 



note one remarkable feature in these POs. One electron, 
by rapidly revolving round the nucleus, screens (or weak- 
ens) the Coulomb field which the other electron feels. 
During this phase, the other electron moves slowly. Two 
electrons alternate in one period between a rapid motion 
near the nucleus and a slow screened motion at a dis- 
tance from the nucleus. Interestingly, the PO, which ex- 
hibits this alternation, is unstable. We found POs with- 
out binary collision through interpolation between two 
POs with binary collision. Some of stable frozen plane- 
tary orbits [17] are among them (see Figs. Qe-g)). 

Finally, we characterize POs as follows. In Figs. [4^a- 
d) , it is clear that one electron revolves round the nucleus 
more than the other. Its number of revolutions up to the 
turn-back point (in the half period) is characterized by 
a half integer or an integer, say n r . For Figs. |4^a-d), 
it is easily inferred that there exists a sequence of POs 
from n r — 1 to n r = oo (i.e., toward ionization). For 
these POs, the Lyapunov exponents and the period are 
numerically calculated (Fig. |5|. It shows a remarkable 
regularity of the period, T = S/2 = an r + 6, where T is 
the period and S is the action. 

Now taking the partial summation over n r in the 
Gutzwiller periodic orbit sum d osc (E) = A p e% Sp with 
a daring approximation, i.e., A p is factored out and is 
replaced by the mean value A, thanks to Poisson sum- 
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FIG. 6: Schematic pictures of binary-collision curves and pe- 
riodic orbits (POs) in -D1/4: (a) Schematic picture of binary- 
collision curves (BCCs) of the first kind and the second kind. 
Solid (resp. dotted) lines represent the BCCs for 2-3 (resp. 1- 
3) collisions. There are infinite number of fundamental blocks. 
The fundamental block is numbered by N from the right hand 
side, (b) POs in the fundamental block N = 1: The triple 
collision orbits in the line a = tt are indicated. The BCCs of 
the first kind are nested. The hexagon (resp. star) represents 
a position of a PO with (resp. without) binary collisions. 
The connected graph, which links the positions of the POs, is 
displayed. 



mation formula, we obtain the Rydberg energy levels, 



En 



c 



with constant C. This means that the regular 
structure found is important for the energy levels of the 
helium atom. The amplitude factor A p may play a role 
of determining the quantum defects for the helium atom 
as observed in the collinear eZe case [4 . To complete this 
program to obtain the energy levels, we have to seek the 
geometrical structure of the whole (unstable) POs. The 
result of the present Letter is a hint to this problem. 

In Fig.[6j a schematic disposition of BCCs and POs are 
illustrated. As shown in Fig.[6ja), there is a fundamental 
block between two BCCs of the second kind. The number 
of fundamental blocks is infinite (Fig. [6^a)). The fractal 
structure of BCCs and the tree structure of POs similar 
to a binary Cayley tree are expected to be topologically 
the same in every block (Fig. |6jb)) reflecting the fractal 



distribution of TCPs on a = n. We show in Fig. [6jb) 
the structure of the rightmost block. We find that there 
exists a PO on each lobe in Fig. [6^b), even on the lobe 
of case (ii). It is interesting to note that though POs are 
far from a = 7r, they keep a kind of one-dimensionality 
along with other orbits in BCCs. 

In conclusion, we have informally proved the existence 
of self-retracing POs with binary collisions, which form 
the backbone of the phase space, using the structure 
of BCCs and have indeed found such POs numerically. 
These findings would serve theoreticians with new stim- 
uli, who would like to carry out the semiclassical treat- 
ment of the helium atom in two dimensions with zero 
angular momentum. 
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